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We present preliminary results from exploring the phase diagram of finite temperature QCD with three 
degenerate flavors and with two light flavors and the mass of the third held approximately at the strange 
quark mass. We use an order a^a^, a** Symanzik improved gauge action and an order agO^, a'^ improved 
staggered quark action. The improved staggered action leads to a dispersion relation with diminished 
lattice artifacts, and hence better thermodynamic properties. It decreases the flavor symmetry breaking 
of staggered quarks substantially, and we estimate that at the transition temperature for an Nt = 8 
to Nt = 10 lattice all pions will be lighter than the lightest kaon. Preliminary results on lattices with 
iVi = 4, 6 and 8 are presented. 

1. INTRODUCTION 

With the Relativistic Heavy Ion Collider (RHIC) now producing data, it has become even more impor- 
tant to understand the phase diagram of QCD at finite temperature, and to determine properties of the 
high temperature quark-gluon-plasma phase with confidence, i.e. with controlled lattice spacing errors. 

It is fairly well established that QCD with two flavors of massless quarks has a second order finite 
temperature, chiral symmetry restoring phase transition. This transition is washed out as soon as the 
quarks become massive. QCD with three flavors of massless quarks has a first order finite temperature, 
chiral symmetry restoring phase transition, which is stable for small quark masses. Not well known is how 
large the quark masses can be before the phase transition turns second order and then into a crossover, 
both for degenerate quarks and especially for the physically relevant case of two light and one heavier 
strange quark. 

In previous studies, the second question is particularly badly answered due to the flavor symmetry 
breaking in Kogut-Susskind quarks, usually used for this purpose: how can one study the influence of 
the strange quark when most of the (non-Goldstone) pions are heavier than the (Goldstone) kaon? 

Adding a few terms to the conventional Kogut-Susskind action, namely three-link, five-link and seven- 
link staples and a third-neighbor coupling, removes all tree- level ©(a^) errors |l|-|. This "Asqtad" 
action shows improved flavor and rotational symmetry ||^,^, and, at least in quenched QCD, good scaling 
properties j^]. This is illustrated in Fig. 0. There, and in other figures below, we plot results in units of 
ri, a scale defined in terms of the static QQ potential by ?'i-FQQstatic(^i) — 1j which leads to ri ^ 0.35 
fm@. 
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Figure 1. The mass of the second Ughtest pion, for fixed fowest pion mass in units of ri, shown as the 
horizontal line, as function of the lattice spacing (left). The "Asqtad" action gives the smallest flavor 
symmetry breaking. The mass of the vector meson, the rho, as function of the lattice spacing (right). 
The "Asqtad" action gives the best scaling. 
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Figure 2. The energy density (left) and pressure (right) of free massless fermions as function of temporal 
lattice size Nf. For free fermions the "Asqtad" action reduces to the Naik action. The "p4" action is the 
improved staggered fermion action preferred by the Bielefeld group ||^ . 



The three-link "Naik" term insures a good dispersion relation and thereby helps decrease the lattice 
artifacts in energy density and pressure at small temporal lattice size Nt, as can be seen in Fig. ^. 

Based on the zero temperature simulations of Ref. Q and the estimate of Tc ~ 150 — 170 MeV |^ we 
deduce that for iVt = 8 — 10 the kaon will be heavier than the heaviest non-Goldstone pion at the finite 
temperature transition. 

We have zero temperature results, in particular, the value of the (bare) strange quark mass, at fixed 
lattice spacing a ^ 0.13 fm and a 0.2 fm. Since we want to keep the physical quark masses approxi- 
mately constant, when we vary the temperature which, at fixed Nt = l/{aT), means varying the gauge 
coupling f3, we interpolated (extrapolated) between the values at the two lattice spacings. 

Given a gauge coupling P and quark mass arrig, for the three- fiavor simulations, or quark masses arus 
and amu,d, for 2 -|- 1 flavors, we determined the tadpole factor uq = (TrC/p/3)^/'' self-consistently in short 
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Figure 3. Real part of the Polyakov line for three flavors with niq « nis (left) and niq ~ 0.6ms (right). 
The data for Nt — 6 have been multiplied by two and the data for A^j = 8 by four. The spatial lattice 
sizes are Ns = 2Nt. 
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Figure 4. (ip'tp) for three flavors with w (left) and « 0.6ms (right). 



zero temperature simulations on lattices. 

To convert lattice scales to physical scales we interpolated a/ri with the form advocated by AUton j|] 

a/n = Cfigl) [1 + C29lf{9l)\ , f{9l) = {bogl)-'^/^''^'^ exp (-^) ■ (1) 

For the (one-loop, tadpole) Symanzik improved gauge action the bare coupling is related to /3 by 
gf) = 10/(3. 6o,i are the universal first two coefficients of the beta-function for three massless flavors, and 
Co, 2 are determined from the measured values of a/ri at a ~ 0.13 fm and a ^ 0.2 fm. 

The Hybrid Molecular Dynamics R- algorithm was used for all simulations, with two noise vectors for 
the 2 + 1 flavor simulations. 
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Figure 5. The ipip susceptibility for three flavors with niq « nis for lattices with Nt = 6 (left) and Nt = 8 
(right). 



2. RESULTS FOR THREE DEGENERATE FLAVORS 

For our first simulation with ruq « we used a linear interpolation of amq{P) in Since the range 
of interpolation/extrapolation is quite large, we subsequently used a linear interpolation of log{amq{P)), 
compatible with the leading behavior of a (with ruq fixed) as function of /?. In the region of the finite 
temperature transition or crossover for lattices with temporal extent iVj = 6 and 8, the difference between 
the two interpolation schemes is small enough that it will not affect our results significantly. 

We show in Fig. || the real part of the Polyakov line and in Fig. ^ the condensate as function of the 
temperature, with the physical scale obtained via ri, for our simulations with three degenerate quarks 
of mass niq « TOs and « 0.6ms. The Polyakov line shows a crossover from confined behavior at low 
temperature to a deconfined behavior at high temperature. The condensate decreases as the temperature 
is increased, but shows no sign of a transition or sharp crossover. 

In Fig. 1^ we show the chiral susceptibility, Xipip = d{'ip^)/dmq from the simulations with niq ~ rng for 
lattices with A^^ = 6 and 8. A peak is seen in both cases. But the peak height does not increase when we 
change the spatial lattice size from Ng = 2Nt to 3Nt , indicating that we are observing a smooth crossover 
from a chirally symmetric phase at high temperature (large /3) to a phase of broken chiral symmetry 
at low temperature. The location of the peak coincides approximately with the steep increase of the 
Polyakov line in Fig. ||. The chiral susceptibility for iriq « 0.6ms looks similar to Fig. H, but in that case 
we have, so far, results only for one spatial size, Ng = 2Nt. 

3. RESULTS FOR 2+1 FLAVORS 

We have also preliminary results from simulations with 2+1 dynamical flavors, with the heavier quark 
kept at the strange quark mass, and the lighter masses at mu,d = O.Gm^ and 0.4ms. The real part of the 
Polyakov line and the light quark condensate from these simulations are shown in Fig. |[ Again we see a 
crossover behavior, particularly for the Polyakov line. For these simulations we do not have the data yet 
to compute the chiral susceptibility. 

4. CONCLUSIONS 

We have made first simulations to explore the finite temperature phase diagram with an improved 
staggered fermion action, the "Asqtad" action, which reduces flavor symmetry breaking so that all pions 
are lighter than the kaon already at larger lattice spacing, and which improves rotational symmetry and 
the dispersion relation, leading to diminished lattice artifacts in energy density and pressure. 

For three flavor simulations down to quark mass m^ ~ 0.6ms, and for 2+1 flavor simulations with the 
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Figure 6. Real part of the Polyakov line (left) and {tpip) for the light quarks (right) for 2+1 flavors. The 
Polyakov line data for Nt = 8 has been multiplied by two. 



light quark masses down to rriu^d = 0.4ms while the heavier quark is kept at the strange quark mass, we 
observed finite temperature crossover behavior, but so far no sign of a genuine phase transition. This 
result is compatible with other recent simulations |]l0| , |rT| which found phase transitions only at quark 
masses lighter than those studied by us so far. The temperature at the crossover, with the scale set by 
ri, Tc ~ 190 - 200 MeV, is a little higher than expected from previous determinations 0,^. This is 
presumably also due to our quark masses still being rather high. 

This work is supported by the US National Science Foundation and Department of Energy and used 
computer resources at Florida State University (SP), NERSC, NPACI, FNAL, and the University of Utah 
(CHPC). 
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